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ABSTRACT 

Sensitivities of nuclear reaction rates to a variation of nuclear properties are studied. 
Target nuclei range from proton- to neutron-dripline for 10 < Z < 83. Reactions con- 
sidered are nucleon- and a-induced reactions mediated by the strong interaction. The 
contribution of reactions proceeding on the target ground state to the total stellar rate 
is also given. General dependences on various input quantities are discussed. Addition- 
ally, sensitivities of laboratory cross sections of nucleon-, a-, and 7-induced reactions 
are shown, allowing to estimate the impact of cross section measurements. Finally, 
recommended procedures to explore and improve reaction rate uncertainties using the 
present sensitivity data are outlined. 

Subject headings: astrophysical reaction rates — nucleosynthesis — compound nucleus 
reaction 



1. Introduction 

Nuclear reaction rates are a central quantity in all investigations making use of reaction net- 
works to study nucleosynthesis and energy generation by nuclear reactions in astrophysical environ- 
ments. The majority of reaction rates appearing in astrophysical plasmas has to be theoretically 
predicted due to a number of reasons. Reactions on highly unstable nuclei, mostly appearing in 
explosive nucleosynthesis, cannot be accessed in the laboratory, yet. Close to and at stability, tiny 
reaction cross sections due to Coulomb barriers often prevent a measurement at astrophysically 
relevant energies. Finally, even if experimental reaction cross sections are available, they are only 
known for nuclei in their ground states. Due to thermal excitation in the plasma, however, a con- 
siderable fraction of nuclei may be in excited states and reactions proceeding on those have to be 
included in the calculation of the reaction rate. 

A frequently encountered question in this context is the one asking for a quantification of 
the uncertainties in the predictions. The knowledge of allowed variations in a given reaction rate 
may be useful to also determine the connected astrophysical uncertainties or to even rule out a 
specific scenario. The question, however, is difficult to address because theoretical uncertainties 
are fundamentally different from the uncertainty quantification generally used as "error bar" in 
measurements. It is possible, however, to identify different factors contributing to the "theoretical 
uncertainty" of a rate and to discuss them separately. 
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The first factor is the choice of the theoretical framework, the model. Obviously, it is a priori 
impossible to exactly determine the correctness of a model without being able to falsify it through 
comparison with data. Even previous knowledge to justify the assumptions can only be used to 
motivate the choice of a model but not strictly prove its validity. This impossibility poses the 
fundamental problem of quantifying uncertainties of the theory approach. It is loosely related to 
systematic errors in experiments due to instruments and external effects, which are also hard to 
estimate. Systematic errors in measurements can usually be removed after identifying the cause. 
For theory, this would mean to give up a model for another one, driven by the failure to describe 
data. This is not just of philosophical importance. As discussed below, an important implication 
of this is the fact that model limitations and applicability can only be estimated and the estimate 
cannot be strictly quantified in a general manner. 

Assuming that a "correct" model is used, the next step in the uncertainty analysis would be 
to check the reliability of the input required to perform calculations with the model. The input is 
comprised of fundamental constants, such as the speed of light or the value of the Planck constant, 
and other quantities, e.g., nuclear masses or other nuclear properties. Part of this input is known 
experimentally and carries well defined statistical errors. These errors can be propagated in the 
model, either through systematic variation within the given limits or by a Monte Carlo analysis, and 
their impact on the final prediction can be quantified. Unfortunately, there are two further classes 
of input, for which this is not possible. One of these classes is experimental input bearing a large 
systematic uncertainty. Specifically for our purposes here, this may include incomplete nuclear level 
schemes, especially from miclcar spectroscopy of unstable nuclides. Finally, some of the required 
input may again be provided by theoretical models. This will be the case for the majority of the 
input required to predict reaction rates for nuclei far off stability. This last input class suffers from 
the same fundamental problem we encounter when trying to quantify uncertainties for our model of 
choice. Therefore it is impossible to assign an " error bar" to this input type and know its impact on 
the final prediction. This cannot be simply remedied by comparing results of different theoretical 
models. This would not be a systematic variation covering the range of possible values. A set of 
models may coincidentally yield very similar values for a given quantity (even if only in the energy 
range relevant for the calculation of a specific rate) but may be equally incorrect. 

The impact of a variation in the input on the final result can be studied independently of the 
actual uncertainties. Of course, this still involves the assumption of having chosen an appropriate 
model to correctly describe the physical situation. Although this will not lead to a full quantification 
of the uncertainties, it provides important information in showing which input quantities have the 
largest contribution to the overall uncertainty in the prediction. This allows to focus on improving 
the uncertainties in the most important quantities and to rule out certain other quantities to have an 
effect on the predicted value. Such an approach also has the advantage that it does not distinguish 
between Tinccrtaintics stemming from experiment and from theory. Having made a choice (or gTiess) 
of uncertainties in the input, they can be simply plugged into the sensitivities of the model to arrive 
at a final error estimate. 
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This work aims at providing tables of the sensitivities of nuclear reaction rate predictions for 
nucleon- and a-induced reactions mediated by the strong interaction. The considered target nuclei 
range from proton- to neutron-dripline for 10 < Z < 83. In addition, sensitivities of laboratory 
reaction cross sections to a variation of the same properties are given. This will aid the conception of 
experiments as well as theoretical investigations to reduce the uncertainties in the most important 
properties. 



2. Definitions 
2.1. Sensitivity 

In order to quantify the impact of a variation of a model quantity q (directly taken froin inpu t 
or derived from it) on the final result Vt, the relative sensitivity ^Sq is defined as ( Rauscheri |20 1 ll ) 



It is a measure of a change by a factor of vq, = i7new/^oid in ^ as the result of a change in the 
quantity q by the factor Vq = gnew/^oid, with ^Sq = Q when no change occurs and ^Sq = 1 when the 
final result changes by the same factor as used in the variation of g, i.e., ^Sq = 1 implies vq, = Vq. 

The sign of the sensitivity ^Sq provides further information on the change and is also given in 
the tables presented in § [H Since both t;^ > and Vq > for the quantities studied in this work, a 
positive sign implies that $7 changes in the same manner as q, i.e., O becomes larger when the value 
of the quantity q is increased. The opposite is true for ^Sq < 0, i.e., 17 decreases with an increase 
of q. 

The above definitions are consistent with the ones used in standard sensitivity analysis when 
realizing that Equation ([1]) is equivalent to 

% = ^ , (2) 

"Jold 

with dfl = rinow - ^old and dq = q^ew - Qold- 

In the following, the final result is either a rate (r) or a reaction cross section (cs). The 
varied quantities are (averaged) widths and the nuclear level density. 



2.2. Reaction rates 

The stellar reaction rate r* is calculated by folding a Maxwell-Boltzmann distribution ^MBiE, T) 
of the relative interaction energies between target A and projectile a with the effective cross sec- 
tion a'^^ which describes the probability that a reaction occurs under the given plasma conditions 
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(|Fowleilll974l : iHolmes et al.lll976l : lRauscheill201ll ') 

* n^n^ ^ ^^^^^^ 

JO 



(3) 



Go(T) 

For simplicity, the reactivity R* (or rate per particle pair) is also called rate in some publications. 
Here, the energy E is measured relative to the ground state (g.s.) of A. The plasma conditions 
are given by the number densities n^, ua of projectiles and targets, and by the temperature T. 
The normalized partition function Gq of the target accounts for the effect of thermally populated 
excited states. 



Go 



(4) 



2 Jo + 1 

with being the spins of ground and excited states, and their excitation energies relative to 
the ground state. (The g.s. is assigned n = and £"0 = 0.) 



^The effective cross section (jFowleii Il974l : IWard &: Fowlerl llQSOl : iHolmes et al.l Il976l : iRauscher 

201 ll ) includes a weighted sum over all possible transition cross sections a^^'^ connecting all initial 
states n in nucleus A at E^ < E with all energetically accessible final states v in nucleus B, i.e., 
with all states with Ey < E + Q, where Q is the reaction Q value, 

£, 



cff 



+ 1 



2 Jo + 1 E 



(5) 



Following [Fowled (|1974| ). the above equation implicitly assumes that cross sections a^^'^ at zero or 
negative energies E — E^ are Zero. 

The effective cross section differs from the definition of the laboratory cross section o"''^^ = 
a^^^ in that it also includes weighted contributions from excited states in the target. This 
accounts for the thermal population of excited states in an astrophysical plasma. Their contribution 
is more important in nuclei with a larger number of low- lying excited states and/or at elevated 
temperature. The g.s. contribution X to the stellar rate. 



provides a measure for the importance of excited states (jRauscher et al.ll201ll ). It can only assume 
values in the range < X < 1, where the laboratory (ground state) cross section is sufficient to 
calculate the stellar rate when X = 1. A smaller X shows increasing contribution of excited state 
transitions to the stellar rate and diminishing importance of the g.s. transitions. For example, 
X = 0.3 implies that 70% of the stellar rate is determined by transitions from excited states in the 
target nucleus and only 30% of the rate stems from reactions on the g.s. The value of X is also a 
measure by how much a theoretical uncertainty can be reduced by a measurement of the laboratory 
cross section (for a further discussion, see 



Stellar rates, but not l aboratory ra t es, obey an iniporta nt reciprocity relation connecting 



forward and reverse reaction (jFowlerlll974l : IWard &: Fowler 



1980) • The stellar reactivity R\j^ of the 
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reaction a + A 
by 



B + b is connected to the reactivity -R^^ of the inverse reaction b + B ^ A + a 



R 



BA 



^AB 



R 



BA 



^AB 



CX 



oc 



rpS/2^-QAB/{kT) 



(7a) 
(7b) 



Equation ([7a|) apphes when projectile a and ejectile b are particles and Equation ([7b|) applies when 
projectile a is a particle and 6 is a photon. The a ,bove relations were deri ved from the well-known 
reciprocity theorem of nuclear reactions (see, e.g. JBlatt Sz Weisskopall952l ) by considering that cj^^ 
includes transitions to all energetically accessible states in the target and residual nucleus. 



Reaction models 



Here we only want to outline the relevant features of the models to be able to di scuss their sen- 
sitivit ie s. For a detailed discussion of reac t ion models used in astrophysics, see, e.g.. iDescouvemont 
(|2003l ): lDescouvemont fc Rauscheil (|2006l ^ : iRausche] (|201lh . 



The majority of nuclear reaction rates were calculated in t he compo und nucleus reactio n 
model, applying the statistical model bvlHauser FeshbachI (119521'). s e e, e.g . , Holmes et al.l (119761): 



Sargoqdl (ll982l'):ICowan. Thielemann. &: TruranI (|l99ll ): lRauscher &: Thielemann 



Wooslev et al.l(ll978l'l : 

(|l998l . l2000l ): lGorielvl (|2003l ') : ICvburt et all (l2010l ). The statistical model assumes a sufficiently high 
nuclear level density (NLD) in the compound nucleus C formed by a reaction a + A^C^B + b, 
with target A, projectile a, residual nucleus B, and ejectile b. In the case of capture reactions, the 
ejectiles are photons. The central assumption in this model is that the compound nucleus is formed 
at an excitation energy where such a large number of resonances of all required spins and parities 
is encountered that they cannot be individually resolved. Instead of decay widths of an isolated 
resonance, the widths of "average" resonances are used. It can be shown that the statistical model 
can be derive d from the single-resonance Breit- Wigner formula by averaging over a large number 
of resonances (jPescouvemont &: Rauscherll2006l ). 



The energy scheme for a compound reaction is shown in Figure [TJ It applies to both the 
case of isolated resonances and the statistical model. The latter would additionally involve a sum 
over many compound resonances Jc, t^c at each Eq whereas there are only few, well separated 
resonances with discrete Eq in the other case. 

The sensitivities of rates and cross sections presented in this work were co mputed in the 
statistical model with version 0.8.3 of the SMARAGD code (jRauscherl |2009| . 



2011 



Nevertheless, 

similar sensitivities would be expected when the used widths were not averaged widths of a large 
number of compound resonances but of a single, narrow resonance dominating the cross section. 
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The rates depend on similar expressions in both cases, a fraction 



HEcJcT^c) = (2Jc + l 



(8) 



The quantities W^, W^, and W^^^ are the total width of a compound level in C with spin Jc at 
excitation energy Ec and including all transitions to levels /i in nucleus A with < Ec — S*] 



pro 



in the entrance channel, the total width of this compound state including all possible transitions 
to levels u in B with E^, < Ec — Se^ in the exit channel, and the combined width W^^ = 
'}2x=A B /c.-Tc f'^Q™ the sum of all emission p rocesses from the compound state to all channels 



A, respectively ([Holmes et al.lll976l : lRauscherll201ll ). The separation energies of the projectile and 



the ejectile in the compound nucleus C are denoted by Spm and Sej, respectively. Assuming narrow 
resonances, the rate of a single Breit-Wigner resonance at excit a,tion energy Eg (this would be an 



energy of E = Eg — Sp^ in the entrance channel) is given by ([Fowler. Caughlan. &: Zimmerman 



19671 : lRauscheijl201lh 



F{Ec,Jc,T^c) 



(9) 



In the Hauser-Feshbach model it is assumed that there is a large number of resonances at any 
compound formation energy Eq and for any Jq and parity -kq- Therefore, the rate is 



Thf oc 



Go(fcr)3/2 



/ V T{E + Sp,o.Jc,T^c)^MB[E,T)dE . 



(10) 



Jc^c 



Although the sum over Jq is unrestricted, only few values will actually contribute, due to quan- 
tum mechanical selection rules and the fact that contributions of large partial waves are strongly 
suppressed by the angular momentum barrier. 

Expressions for the laboratory rates are simply obtained by dropping the factor I/Gq in Equa- 
tions ([9|) and (jlOp and by using only one transition in t he determination of the appea ring in 
the numerator of Equation ([8|), the one to the g.s. of A ([Holmes et al.lll97d : iRauscheillioill ). 



The importance of direct reactions - and specifically of di rect neutron capture - in astrophysical 



envir o nments has also been discussed i n literature (see, e.g., iMathews et al.lll983l : iRauscher et al 



19981 : iGorielv. Hilaire. &: Chamell l2009l : iRauscherl l201ll . and references therein) . These are not 
included in the present sensitivity study. In Figure [J they would be shown as positive-energy 
transitions directly between thermally populated states in the nuclei A and without the formation 
of an intermediate compound state in C. The energy difference between initial and final states 
would be emitted in a single photon. Their sensitivities, however, are trivial to estimate because 
only single transitions are involved instead of summed widths. Thus, the rate of a direct reaction 
will linearly depend on the spectroscopic factors used, with a sensitivity s = 1. They will also be 
more strongly sensitive to a change in the properties of the single final state (excitation energy, 
spin, parity) than compound reactions ([Rauscher et al.lll998l ). 
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Direct reactions at high energies ~ (d,p) and (d,n) reactions at several tens to hundreds of MeV, 
for instance - have also been suggested to be used in nuclear spectroscopy studies of unstable nuclei 
at radioactive ion beam facilities. Cross sections of reactions at such high energies do not appear 
in the calculation of astrophysical reaction rates but the extracted information on low-lying excited 
states and spectroscopic factors is useful in the calculation of the widths W implicitly appearing 
in Equations dH) and ([lO]) (see also § [42]l . 



3.1. Model applicability 

As pointed out above, the statistical model relies on the assumption that a large number of 
narrow compound resonances is present at the formation energy of the compound nucleus. Whether 
this is valid needs to be ascertained, in principle, for each reaction and incident projectile energy 
separately. As this would involve a comparison of the actual number and widths of resonances to 
the averaged widths predicted by the model, this is not strictly possible within the model itself. 
This is the first type of fundamental uncertainty described in § [TJ It has to be mentioned, however, 
that there is a difference in the applicability of the Hauser-Feshbach model to reaction cross sections 
and to astrophysical reaction rates, as these involve different energy distributions of the projectiles 
and therefore average in a different manner. In an experimental determination of a cross section 
the energy distribution of the projectiles is given by the energy distribution within the original 
particle beam and the energy straggling within the target. Considerable efforts are usually made 
to confine the resulting energies to a narrow range, in order to limit the energy uncertainty in the 
measured cross section. Therefore, resonances with a width comparable to the obtained energy 
resolution will be noticeable in the excitation function. In the calculation of a reaction rate, on 
the other hand, the cross section is averaged over a Maxwell-Boltzmann energy distribution (see 
Equation [3]) which is much wider than the energy resolution of an experiment. This explains the 
fact that the reaction rate is rather "forgiving" to deviations around the "true" cross section value, 
provided the deviations cancel within the integration. In consequence, the Hauser-Feshbach model 
is more widely applicable for reaction rates than for cross sections. It will even be applicable in the 
presence of narrow resonances within the astrophysically important energy range as long as their 
average contribution within this range is correctly accounted for. 



Applicability maps for neutron-, proton-, and a-induced reactions were shown in lRauscher. Thielemann. &: Krs 



(|l997l ). specifying the lowest stellar temperature at which the statistical model is estimated to be 



still suitable to predict the stellar rate. (It should be noted that these were estim ates based on a pre- 



dicted NLD.) The temperature limits were also explicitly included in the tables of lRauscher &: Thielemann 



(I2OOC1 ). The rates for the majority of reactions with Z > 10 can be described within the statistical 
model. The exceptions are rates for which the dominant contributions to the reaction rate integral 
stem from compound excitation energies with low NLD. For neutron- induced reactions this can 
be the case even close to and at stability for closed shell nuclei. The Gamow windows are shifted 
to higher compound excitation in reactions for which the energy dependence is determined by a 
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charged-particle width. Therefore these are not so sensitive to closed shells at stability. Reactions 
on heavy nuclei, however, may involve a-induced reactions with negative Q value. This may move 
the Gamow window back into a region of low NLD. Moreover, proton- and neutron-induced reac- 
tions close to the proton- and neutron dripline, respectively, also show low or negative Q values 
and also for these cases the statistical model cross sections may be suited only to calculate rates at 
high temperature. 

The astrophysically r elevant energy windows as function of temperature have recently been 
revised bv iRauscherl (l20ld'). This impacts the derived applicability limits only in a limited manner. 



as explained in § 5.3.2 of iRauscheij (|201ll ). 



4. Results and Discussion 

4.1. Sensitivity tables for width variations 

4.1.1. Rates and ground state contributions 

Table [J contains the sensitivities (as defined in Equation [T|) for astrophysical reaction rates 
when systematically and separately varying the neutron-, proton-, a-, and 7- widths by a factor 
Vq = 2. These widths are always included in VF^" of Equation (llOp and, depending on the reaction 
studied, can also appear as W^^^ or W^^. The g.s. contribution X to the stellar rate is also given. 
The table is organized as follows: The first column contains the abbreviated element name of the 
target nucleus A and the next two columns give its charge and mass number Mj^. The fourth 
and fifth columns identify the reaction by specifying projectile a and ejectile 6, respectively. Roman 
letters are used in the identification of projectile and ejectile: n (neutron), p (proton), a (a), g (7). 
This is followed by 24 blocks for the 24 temperatures 0.1, 0.15, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 
1.0, 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, 5.0, 6.0, 7.0, 8.0, 9.0, and 10.0 GK. Each block contains five 
columns showing the sensitivity ^Sq to a variation of the 7-, neutron-, proton-, and o-width {^Sg, 
'^Sn, '^Sp, ''sa, in that order) and the g.s. contribution X (Equation [6]) at the given temperature. The 
table is sorted by Z, A, and reaction type. 

Except for captures, only reactions with positive reaction Q value are shown. Due to the 
reciprocity of stellar rates (Equations \7a\ \7h\i , the sensitivities are the same for forward and reverse 
reaction. The g.s. contributions, on the other hand, can be quite different. For the majority 
of rates, X is smaller for the reaction direction with negative Q value. Only comparatively few 
exceptions are known to this rule. They are due to a Coulomb suppression eff^ect of excited state 



contributions which is discussed in iKiss et al.l (|2008l ): iRauscher et al.l (120091 ). Because it may be 



of interest for experimental nuclear physics, the g.s. contributions X for reactions with negative Q 
value are given in Tabled A similar table (Table [3|) is given for photodisintegration rates. This is 
mainly to illustrate that especially photodisintegration rates have much smaller g.s. contributions 
than the respective capture rates. In consequence, the determination of a photodisintegration 
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cross section on a target nucleus in the g.s. does not provide much information on the stellar 
rate. Moreover, it exhibits a different sensitivity to the input as can be seen from the tables 
provided in ^ I4.1.2I To provide a quick overview, Figures [2H7] show X at 1.5 GK for captures and 
photodisintegrations involving neutrons, protons, and a-particles. The g.s. contributions will be 
smaller at higher temp eratures. Value s of X for (n,7) close to stability at s-process temperatures 



have been discussed in Rauscher et al. (2011 



The general behavior of the sensitivities can be understood by recalling a few simple facts. 
First, it has to be noted that in astrophysical applications forward and backward rates are in equi- 
librium at temperatures above 5-6 GK. Above this temperature, individual rates are not important 
and the attained equilibrium abundances are governed by the reciprocity relations in Equations 
(j7ap and (|7bp . The temperature dependence of the sensitivity of a reaction rate to a given width is 
comparatively weak below 5 GK. This is because the relevant energy range covered in the integra- 
tion for the rate (Equation [3|) is comparatively narrow as it depends on the width of the Gamow 
peak at the given temperature and the effective cross section is integrated over this range. It is 
even narrower for neutron channels where the relevant energy range is given by the width of the 
Maxwell-Boltzmann distribution. 

As for the dependence of the rates on certain widths, a few cases can be distinguished regarding 
the ratio of widths in Equation ([8]) : 



1. The fraction foremost depends on the smallest width appearing in the numerator when 
the value of the denominator is dominated by the larger width in the numerator. This is 
frequently the case as the denominator contains the sum of all channels including the ones 
appearing in the numerator. In this case the denominator cancels with the larger width and 
the smaller width remains. Its energy dependence then determines the energy dependence of 
the cross section. The rate and cross section will then be completely insensitive to changes in 
the larger width (as long as it does not become comparable or smaller than the other width 
in the numerator) and any change in the smaller width will transfer fully to the cross section 
and rate. 

2. When the denominator is not dominated by either of the two widths appearing in the nu- 
merator, it is irrelevant whether they are equal or of different magnitude. The cross section 
and rate will be sensitive to a change in any of the widths, either one in the numerator and a 
third width dominating the denominator. A change in the denominator, however, will affect 
the rate in inverse proportionality, i.e., an increase in a width in the denominator will lead to 
a decrease in the cross section and rate. 

3. The most complicated case appears when the two widths in the numerator are comparable 
and both are strongly contributing to the total width in the denominator. The cross section 
and rate will then be sensitive only to changes in these two widths but the level of sensitivity 
depends on the size of the variation factor. The sensitivity will be low with a small variation 
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factor because any change in the numerator is cancelled by a similar change in the denomi- 
nator. Using a large variation factor, on the other hand, may transform this case to case 1 
above, when the varied width becomes much larger or smaller than the other ones. 



Case 1, of course, is well known in resonance analyses performed in experimental nuclear 
physics. A closer inspection of the sensitivities, however, shows that a rule-of-thumb often applied 
in nuclear reaction analysis does not always work for astrophysical rates. The rule states that the 
7" width is smaller than the particle widths and therefore the energy dependence and sensitivity of 
a capture rate will be solely governed by the 7-width, according to case 1 above. Astrophysical 
charged particle capture on intermediate and heavy nuclei (with > 20), however, shows a 
different behavior. This is due to the fact that the interaction energies of the charged particles 
are close to or below the Coulomb barrier and consequentially the widths are suppressed. They 
become comparable or smaller than the 7-width and thus determine the sensitivity behavior. This 
implies that 7-widths play a minor role in certain astrophysical charged-particle capture reactions, 
contrary to the case of standard nuclear reactions at higher energy. They still play a role in neutron 
captures, though, unless the relevant energy window is close to the neutron threshold. Close to the 
threshold, the neutron width may become smaller than the 7-width again. 

Figures [8] and [9] illustrate the above by showing the sensitivities on different width variations 
of the stellar rates for "'^'^^Sm(a,7)"^^^Gd and "^^^Sn(n,7)"'^^^Sn, respectively. While the a capture 
reaction is only sensitive to the a-width which remains smaller than the 7-width across all astro- 
physically relevant temperatures and energies, the neutron capture reaction is indeed governed by 
the 7-width. It should be noted that the neutron capture Q-value of the doubly magic nucleus 
^^^Sn is small and therefore the NLD at the compound formation energy may be too small to 
allow application of the statistical Hauser-Feshbach model (see § 13. ip . Single resonances and direct 
capture may provide important contributions to the reaction rate. 

As another example, the sensitivities of the stellar ^^Ru(p,7)^^Rh rate are shown in Figure 
[TOl The rate is mostly sensitive to the proton width over all of the relevant temperature range. 
Towards higher temperatures, however, the sensitivity to the proton width is slowly decreasing 
while the one to the 7 width is increasing. At 5 GK both sensitivities already are quite similar: the 
sensitivity to the proton width is only ^Sp = 0.5 while the one to the 7 width is already = 0.3. 
Uncertainties in other widths do not play a role. 

A more general view is provided in Figures [HHISl showing the values of sensitivities across the 
nuclear chart for a given reaction, width, and temperature. Sensitivities for proton capture rates 
at T = 1.5 GK are shown in Figures [TT] and [T2l As seen in Figure [TTl the (p,7) rates are maximally 
sensitive to changes in the proton width all across the nuclear chart except for the lighter range of 



nuclei at th e proton dripline. These have very low reaction Q-valu es and the statistical model is not 
applicable ]R„.„.ch„ Thi.l.n,an.. K^.tB E;;:;;^ ^ . Their stellar rate sensitivity is 



dominated by the sensitivity to the 7-width, as depicted in Figure [T2l For proton-rich nuclei with 
small to moderate capture Q-values (up to a few MeV), the sensitivities to proton or 7-widths are 
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complementary, as a comparison of the two figures reveals. For proton captures on neutron-rich 
nuclei, the rates are sensitive to both widths because the total width in the denominator of J- in 
Equation ([8]) is dominated by the neutron width. This is case 2 discussed above. Therefore they 
are also sensitive to, and change inversely proportionally to, a variation of the neutron widths in 
that region of the nuclear chart. 

A somewhat similar behavior is found for (0,7) rates. Due to the higher Coulomb barrier, and 
thus overall smaller width, the rates are sensitive to variations of the a-width all over the nuclear 
chart. Therefore we only show the sensitivities to changes in the 7-width in Figure [T3l Similar to 
proton captures, a-capture rates on neutron-rich nuclei are sensitive to a-, 7-, and neutron widths, 
with the neutron width dominating the denominator and thus leading to an inverse proportional 
dependence on variations of the neutron width. 

The sensitivities of neutron capture reactions at 0.3 GK (close to typical s-process temper- 
atures) to variations of the neutron- and 7-widths are shown in Figures [U] and [T5l respectively. 
Again, the sensitivities are basically complementary, with decreasing sensitivities to the neutron 
width with decreasing Q-value towards the neutron-rich side of stability. The exception are (n,7) 
rates involving target nuclei close to the proton dripline. These are sensitive to both, neutron and 
7-widths, in the same manner and inversely sensitive to the proton width which determines the 
denominator of as can be seen in Figure [TBI 



4-1.2. Cross sections 

Although the laboratory cross section, i.e. the cross section of a reaction proceeding on a target 
nucleus in the g.s. as defined in § 12.21 cannot be used to calculate the stellar reaction rate unless 
X = 1, experiments may be performed to better constrain input required for the prediction of 
the rate. This is only possible when the measured cross section is indeed sensitive to the nuclear 
property in question. For this purpose, tables of laboratory cross section sensitivities to the same 
variation of the same quantities as for the rates shown in § 14.1.11 are provided here. Table H] is 
structured similarly to Table [T] but it contains the sensitivities to the four widths at a varying 
number of projectile energies and no distinction is made regarding the sign of the reaction Q value. 
Also, the g.s. contribution X is not given, as it does not apply here. Table H] is organized as follows: 
The first column contains the abbreviated element name of the target nucleus A and the next 
two columns give its charge Za and mass number Ma- The fourth and fifth columns identify the 
reaction by specifying projectile a and ejectile b, respectively. The sixth column gives the number 
of energies nEn for which values are present. This is followed by 50 blocks of 5 columns each. 
In each block the first entry is the center of mass projectile energy E in MeV and the next four 
entries are the sensitivities to a variation of the 7-, neutron-, proton-, and a- width C^^Sg, ^^s^, 
''^^Sp, ''^^Sa, in that order). Empty blocks, i.e. blocks beyond n-En, contain an energy of —1.0 and all 
sensitivities are zero. Empty blocks appear when the number of valid energy entries is smaller than 
50. Valid energies are those below 14 MeV for which at least one of the sensitivities is non-zero. If 
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there is a sequence of energies with only zero sensitivity, e.g. for endothermic reactions, only the 
entries for the highest of these energies is shown to indicate that all sensitivities are zero below 
this energy. Reactions with fewer than five valid energies are omitted. The energy cut-off at 14 
MeV still encompasses all astrophysically relevant energies and at the same time ensures that the 
reaction model used is sufficient. 

The cases discussed in § 14.1.11 regarding the dependence of rates on the widths appearing in 
the factor T also apply to cross sections. The energy dependence of the sensitivities, however, are 
stronger than the dependence on the temperature in the rates. The energies shown in Table H] 
cover a larger range than required for the integration of the reaction rates. Moreover, the energy 
dependence in the cross sections is not diluted by an integration as for the rates. In particular, large 
changes in the sensitivities appear at energies where a new reaction channel opens. This underlines 
the fact that it is necessary to experimentally determine the cross sections as close as possible to 
the astrophysically relevant energy range because the cross sections are most likely depending on 
different widths at higher energies. Again, this can be best seen when examining charged-particle 
capture. At high energies, the energy dependence of the cross section and its uncertainty is governed 
by the 7-width whereas in the astrophysical energy range it is determined by the charged-particle 
width. 

For illustration, Figure [TT] shows the sensitivities of the laboratory reaction cross section of 
^^Ru( p,7)^^Rh as functi on of cm. energy. This reaction appears in p-process studies (jArnould &: Goriely 



2003 



Zhong et al.ll2010l ) and the astrophysically relevant energy range is 1.5 < E < 3.85 MeV (for 
plasma temperatures 2 < T < 3 GK). It can be clearly seen that the impact of each width type 
drastically changes from lower to higher energy and thus also the relevance of each type of width 
strongly depends on energy. At low energy, the cross section is only sensitive to the proton width 
and all uncertainties therein fully affect its prediction. At the high end of the shown energy range, 
the impact of the proton width is negligible and the by far highest sensitivity is to the 7-width. 
Uncertainties in the a-width are irrelevant for this reaction over the whole energy range, whereas 
the importance of the neutron width quickly increases above the (p,n) threshold with increasing 
energy. Since the neutron width appears in the denominator of J-" in Equation ([8]), the size of the 
cross section varies inverse proportionally to the size of the neutron width. If a comparison between 
experimental data and theory finds a discrepancy above, say, 7 MeV, this is still inconsequential 
for the astrophysical rate. On the other hand, if theory reproduced the cross sections above 7 MeV 
well, the astrophysical rate would still be unconstrained because it is sensitive to a different nuclear 
property. 

Not all laboratory cross sections have strongly varying sensitivities with ener gv. Figure [181 



shows the situation fo r ^^Ge(n,p)^^Ga, important in i^p-process nucleosynthesis ( Frohlich et al 



20061 ; iPruet et al.ll2006l ). It is exclusively sensitive to the neutron width across the shown range of 
energies. 

As another example, the cross section of ^^'^Hf(n,7)^^^Hf seems to be mostly sensitive to the 
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7-width across the considered energy range. A closer inspection, however, shows that there is a 
quite similar sen sitivity to both 7- an d neutron width at typical energies of — 100 keV encountered 
in the s-process (jGallino et al.l Il998l ) . This is shown in Figure [T9j The laboratory cross section of 
the reverse reaction ^^^Hf(7,n)^^*'Hf, on the other hand, exhibits only sensitivity to the 7-width, 
as can be seen in Figure [20j Although there is reciprocity between s ingle transitions according to 
the reciprocity theorem of nuclear reactions (iBlatt &: Weisskopall952l ). laboratory cross sections do 
not obey reciprocity because different transitions are included in the forward and reverse reaction, 
as can be easily inferred from Figure [T] and also from the definition of the laboratory cross section 
provided in 



4.2. Dependence of widths on input quantities 



The tables in § 14.11 allow to judge the importance of a width for the stellar reaction rate or 
the laboratory cross section. It is then important to consider which nuclear pro perties det e rmine 



a width and its uncertainty. We do not want to repeat the detailed discussion in iRauscheii (j201ll ) 
here but rather give a short summary of what enters the determination of a width. It equally 
applies to the widths of individual resonances to be used with Equation ([9]) and to the averaged 
widths appearing in Equation ()10p . 



The transition widths in Equation 



Holmes et al. 


1976: 


Cowan. Thielemann. & Truran 


1991 



are actually derived by accounting for all ener- 
-7, Jc,vrc) to levels in channel A (see, e.g., 
, and references therein). 



^last 



X 



, n=0 




Ec — S-pr. 



Wil^C,Jc,nc)^{E\j\n^)pX^^X^ jX^ ^A) 



(11) 



Optical potentials are required to compute the particle transition probabilities by solving the (time- 
independent) Schrodinger equation. Photon strength functions, on the other hand, determine the 
strengths of the 7-transitions. In addition, the calculated widths depend on the discrete level 
scheme included up to the n^^^^-th level at energy E^^^^. They may also depend on the NLD 

lA 

last 



when transitions to states above E^\^^ are possible (see Figure d]) 



In principle, the NLD enters in two ways. Firstly, the NLD at the compound formation energy 
determines which reaction model should be used (see § 13. ip . Once the model has been chosen, 
however, it does n o t ente r the calculation anymore (unless using a modified, non-standard model. 



see, e.g.. 



Rauscheii (|201ll )). As discussed in ^ 13.11 the choice cannot be made from within the model 
and therefore comprises a type 1 uncertainty which cannot be estimated systematically as explained 
in § [TJ Secondly, the NLD below the compound formation energy may be required to calculate the 
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widths W appearing in Equations ([9]) and (jlOp when transitions to unresolved or unknown excited 
states are important, as shown in Equation (jlip . Because NLD studies often take their motivation 
from an apphcation to astrophysics, the sensitivities of reaction rates and cross sections to variation 
of the NLD is studied in detail in the following § 14.31 

Nuclear masses are important insofar as mass differences determine Q values and separation 
energies and thus also the relative energies of the possible transitions, as shown in Figured] This 
then impacts the calculation of the transition widths. Therefore, uncertainties in masses become 
relevant when they lead to sizeable changes in separation energies, exceeding a few tens of percent. 
Close to stability the mass differences are known with such high accuracy that essentially no 
uncertainty is introduced into the rates from masses. This holds as long as separation energies 
are large. Approaching the driplines, separation energies become smaller and uncertainties in 
masses larger. Nevertheless, whenever experimental information is available, the uncertainties 
are negligible. For most nuclei along the dripline, only theoretical mass predictions are available 
and different mass models yield largely different values. It should be noted, however, that mass 
differences are still thought to be predicted more accurately than masses of single nuclei. 

Sensitivities of rates to a variation of masses are not shown here for three reasons, the first 
reason being the mentioned fact that the separation energies for the majority of reactions in the 
nuclear chart are known with sufficient accuracy. The second reason is that masses cannot be varied 
independently as they simultaneously affect the separation energies appearing in several reactions. 
The third reason is that for reactions with low separation energies and large mass uncertainties the 
ratio between forward and reverse rate is more affected by the mass uncertainties than the forward 
rate itself. This can be seen in Equations ((Taj) and (|7bp . where the exponential dependence on the 
reaction Q-value is explicitly shown. In most cases, the change in the balance between forward and 
reverse rate will affect astrophysical applications more than the pure change in the reactivity. This 
is only different very close to the driplines, when a change in the masses leads to a largely different 
range of relative channel energies and consequently of accessible levels. Then a recalculation of the 
reactivity is called for. 

4.3. Sensitivity tables for NLD variations 

4.3.1. Rates 

Table [5] shows the sensitivities of the astrophysical reaction rates to variations in the NLD by 
a factor Vq = 2. The table is constructed similarly to Table [U but with several additional columns. 
The first column contains the abbreviated element name of the target nucleus A and the next 
two columns give its charge and mass number Ma- The fourth and fifth columns identify the 
reaction by specifying projectile a and ejectile b, respectively. The sixth to tenth column gives the 
excitation energy -Eiast of the highest-lying, discrete excited level included in the calculation for the 
four reaction channels defined by emission of 7 (g), neutron (n), proton (p), or a-particle (a). This 
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is followed by 24 blocks for the same 24 temperatures as in Table [TJ Each block contains eight 
columns showing the sensitivity ^Sq to a variation of the NLD in the 7-, neutron-, proton-, and 
a-channel. The first set of four columns in a block show the sensitivities (^Sg, j^s^, i)Sp, ^Sa, in that 
order) when including experimentally known excited states in each channel up to the excitation 
energy Ef^^, -^lasf -^last' -^lasf respectively. The second set of four columns in a block give 
sensitivities (^Sg, ^5^, ^Sp, ^s^, in that order) calculated without excited states, i.e., using the 
NLD immediately above the ground state. 

It should be noted that the shown sensitivities actually arise from a combination of two sensi- 
tivities, although they were directly inferred from the variation of the NLD and its impact on the 
final rate, following Equation ([1]). The sensitivity of a width to a change in the NLD is folded with 
the sensitivity of the rate to a change in the width. The latter sensitivity is discussed in § 14.1. H and 
shown separately in Table [H Obviously, when a rate is not sensitive to a change in a width, it will 
also not be modified when this width is changed by a variation of the NLD. As discussed in 14.1.1^ 
it should also be noted that a relevant width can appear in the numerator or denominator of the 
ratio J- in Equation If it is in the denominator, the rate will change inversely proportional to 
the change in the width, i.e., it will become smaller when the width is increased. As before, this is 
identified by a negative sign of the sensitivity given in the tables here. 

The standard approach in the calculation of reaction rates for astrophysics in a Hauser- 
Feshbach statistical model is to include a number o f discrete excite d states, when experimentally 



known. The discrete levels used are a selection from lENSDFl (j201Cl ). including levels with well de- 
termined spins and parities and up to the energy at which the level scheme can still be considered 
complete. A predicted NLD is only used above the energy of the last included experimental level 
Eiast- Figure [2T] shows -Eiast for all considered nuclei. This reduces some of the rate uncertainty 
but obviously only for nuclei with accurately known level schemes. Moreover, this procedure can 
only reduce the uncertainty in the channels where transitions to the included states dominate the 
rate. This occurs mostly in the particle channels but is not the case for the 7-channel of a reaction 
du e to the h i gher Q- value and thus larger range of possible relative energies. It has been shown 



m 



Rauscherl (120081 ) that the dominant 7-transitions in astrophysical capture lead to levels 2 — 4 
MeV below the formation energy of the compound nucleus. Even close to stability, complete level 
schemes are not known at this excitation energy. In captures with low Q-value, on the other hand, 
transitions to the lowest excited states are favored and the knowledge of a few levels above the 
compound g.s. will be sufficient. To see the impact of the included levels, the sensitivities ^Sg and 
pSg, with and without experimental levels, respectively, can be compared for a specific reaction 
rate at a given temperature. The difference will be larger, the larger E + Q — ii^iast is in the reaction 
channel, for channels with -Blast > and Q being the reaction Q-value applicable to the channel. 
The energy E mostly contributing to the reaction rate at a given temperature is the location of t he 



maximum contribution within the relevant energy window and can be found in iRauscherl ( 2010l ). 



As an example. Figure [22] shows the stellar rate sensitivities of ^^Ru(p,7)^''Rh to a variation of 
the NLD in the proton- and 7-channel, with and without inclusion of discrete excited states. This 
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is to be compared to Figure [T0\ showing the sensitivity of this rate to width variations. Firstly, 
it can be seen that the inclusion of excited states does not change much in the NLD sensitivity 
here. Secondly, the rate is not sensitive to the NLD variation in the proton channel although 
it shows the highest sensitivity to the proton width. The sensitivity to a NLD variation in the 
7-channel, on the other hand, is the same as the overall sensitivity to a variation of the 7-width. 
This can be understood by examining which transitions are important. As pointed out above, the 
most important 7-transitions from the compound state end at levels with high excitation energy 
and thus the NLD directly determines the size of the 7-width. In the proton channel, the most 
important transitions go to low-lying states or, as in this case, to the g.s. This fact is confirmed by 
a value of the g.s. contribution to the stellar rate X ~ 1. 

Similar considerations apply for the rates of "'^^^Sm(a,7)"'^^®Gd and "'^'^^Sn(n,7)"'^'^^Sn. Their 
sensitivities to width variations were discussed in § 14.1.11 and shown in Figures [8] and [H Both rates 
are found to be insensitive to variations of the NLD. The rate of ^^^Sm(a,7)^^^Gd is only sensitive 
to the a-width and has a strong g.s. contribution to the stellar rate. This implies that a transitions 
to the g.s. of ^^^Sm are dominating which makes the width insensitive to a variation of the NLD at 
higher excitation energy. Although the rate of ^^^Sn(n,7)^^^Sn is almost exclusively sensitive to the 
7-width, it shows also a very low sensitivity to any NLD variation. The dominating contribution 
to the 7-width stems from the g.s. 7-transition due to the low reaction Q-value. 

In summary, a variation of the NLD only affects rates when the NLD is changed in a channel 
to which the rate is sensitive (see § I4.1.ip and transitions to higher lying states are important. 
Therefore the rates are insensitive to NLD variations in the entrance channel of charged particle 
captures throughout the nuclear chart and only weakly sensitive in neutron captures. They are 
sensitive to a NLD variation in the 7 exit channel when they are sensitive to a variation of this 
7-width at all. However, in these cases charged particle captures are also inverse proportionally 
sensitive to a variation of the neutron width, which enters the total width in the denominator of 
Equation ([8]), and therefore also to a NLD variation in this channel. Reactions with particles in all 
channels exhibit only a very weak sensitivity to NLD variations and this only close to stability for 
Za ^ 65, if at all. 

4-3.2. Cross sections 

Similarly as for the cross section sensitivities to width variations, shown in § 14.1.21 Table [6] 
provides the sensitivities of laboratory cross sections when systematically varying the NLD. As for 
the rate sensitivities in Table [5l these are shown with and without the inclusion of discrete excited 
states. These excited states were the same as those used in Tableland therefore the values for £^iast 
remain the same and are not specified anymore. The structure of Table [6] is close to that of Table 
m The first column contains the abbreviated element name of the target nucleus A and the next 
two columns give its charge Za and mass number Ma- The fourth and fifth columns identify the 
reaction by specifying projectile a and ejectile 6, respectively. The sixth column gives the number 



-17- 



of energies nEn for which values are present. This is fohowed by 50 blocks of 9 columns each. In 
each block the first entry is the center of mass projectile energy E in MeV and the remaining set 
of eight columns contains the sensitivities to a variation of the NLD in the 7-, neutron-, proton-, 
and a-channel. The first four columns in this set show the sensitivities (j^Sg, gsn, ^'^p) D^a, in that 
order) when including experimentally known excited states in each channel. The next four columns 
give sensitivities C^s^, gs^, gsp, gs^, in that order) calculated without excited states, i.e., using 
the NLD immediately above the ground state. As for Table HI empty blocks contain an energy of 
—1.0 and all sensitivities are zero. Photodisintegration and other reactions with negative reaction 
Q value are included. 

The sensitivites of laboratory cross sections to NLD variations display similar behavior as 
outlined above for the rates. In addition, it has to be remembered that the entrance width of a 
reaction appearing in the numerator of Equation ([8]) only contains the g.s. transition and therefore 
is not sensitive to the NLD at all. The total width in the denominator, however, contains the full 
width of the entrance channel, including possible contributions from the NLD as shown in Equation 

As a specific example. Figure [23] shows the sensitivities of the laboratory cross section of 
^^Ru(p,7)^^Rh to a variation of the NLD in the various channels, with and without inclusion of 
experimental excited states. The variation in the a+^^Tc channel is not shown because the cross 
section is insensitive to it. The inclusion of excited states mostly affects the neutron channel above 
the threshold and it also moderately affects the proton channel above 5 MeV. When using the 
discrete excited states, there is no sensitivity to a change in the NLD in both channels. This 
implies that a sufficient range of excitation energies is covered by the experimental spectroscopic 
information. The sensitivity in the 7 channel is only weakly modified. Comparing to Figure [T7] it 
is apparent that transitions to excited states do not play a significant role in the proton channel 
below 5 MeV. The NLD variation does not result in a change of the cross section at those energies 
although the cross section would be very sensitive to changes of the proton width. Interestingly, 
above 5 MeV the cross section variation goes the opposite way as the NLD variation. This is due to 
the fact that the entrance width in the numerator of Equation ([8]) is independent of any NLD but 
the full entrance width appearing in the denominator depends on it to a certain fraction. Obviously, 
such a behavior is never found in stellar rates. Finally, the sensitivity in the 7-channel to the NLD 
variation is similar to the one of the cross section to the variation of the 7-width, showing that 
7-transitions to higher lying excited states in ^^Rh are dominating this width. 

As for the further examples given in the preceding sections, neither stellar rate nor laboratory 
cross section of ^^Ge(n,p)^^Ga are sensitive to NLD variations. The laboratory cross section of 
^^^Sm(a,7)^^^Gd above 12.3 MeV is weakly sensitive to the NLD in the 7+^^^Gd channel but 
insensitive to a NLD variation in the a+^^^Sm channel although this is dominating the cross 
section (and rate) behavior (Figure [M]) • Like before, the reaction ^^'^Hf(n,7)^^^Hf exhibits a more 
complicated behavior at low energy. Comparing Figures [25] and [19] it is found that a NLD variation 
in the 7-channel leads to the same variation as the variation of the 7-width above 40 keV. Below this 
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energy, the sensitivity to the NLD is weaker, indicating that 7-transitions to discrete states are of 
some relevance. Interestingly, the dependence on the NLD in the neutron channel is quite different 
from the sensitivity to the neutron width. While the cross section below 200 keV is insensitive 
to the NLD variation in the neutron channel, implying a domination of transitions to the g.s., it 
changes opposite to the NLD variation at higher energy. This can again be explained by the fact 
that the full neutron width is only included in the denominator of T when dealing with a laboratory 
cross section. 

4.4. Discussion of application 

The tables presented in the preceding sections contain the necessary information to assess 
uncertainties in astrophysical reaction rates. The data presented are of importance also for the 
conception of theoretical and experimental studies to reduce said uncertainties. Table [7] summarizes 
the conclusions from § U] on what input quantities impact the widths. 

As explained in the introduction, the error in the choice of models for both the reaction rate 
calculation and the input to the calculation cannot be quantified. Therefore it also depends on the 
type of input whether a systematic variation can be performed and what inherent uncertainty is 
assigned to the input. A further discussion goes beyond the scope of this work. Nevertheless, only 
a separation of the treatment of input uncertainties and rate sensitivities allows further insight in 
the dependences and their relevance. Although uncertainties in the input have to be determined 
separately, the sensitivities provided here are indispensable in determining the final impact of such 
uncertainties. Knowing that a certain input quantity q carries a given uncertainty factor Ug, with 
the actual value Vq of the quantity being in the range Vq/uq < Vq < UqVq, Table [7| shows which width 
c in the calculation is affected by it. This helps to estimate the uncertainty factor Uc of c which is 
often close to Uq (except for masses). The sensitivities to a variation of that width can be found 
in Table [H The resulting uncertainty factor U in the astrophysical rate then is U = 'scUc- In the 
case of a known systematic uncertainty factor in the NLD in channel c. Table [5] can directly be 
used to obtain J^ScuJ?. 

Practical applications will usually first determine the range of temperatures relevant to nuclear 
burning and a number of important reactions in a specific astrophysical process. Table [T] is then 
used to understand on which widths the astrophysical rate actually depends. The uncertainty in 
the relevant input and how it propagates to the width has to be estimated subsequently. Finally 
the resulting uncertainty in the rate is again computed with the sensitivity found in Table [TJ 

Experimentalists interested in improving a specific rate, i.e., reducing its uncertainty, should 
also first note the plasma temperature range in which this reaction contributes to a given astrophys- 
ical process. The g.s. contributions X in this temperature range, taken from Tabled! allow to judge 
whether a measurement of the laboratory cross section actually can constrain the astrophysical re- 
action rate. This is only the case for X being close to unity. In fact, the value of X directly shows 
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the maximally possible reduction of a rate uncertainty factor U > 1. Assuming that the theoretical 
uncert ainty factors of ground and excited state contributions are similar, U = Uf^ = U^^'^, equation 



(18) of lRauscher et al.l (j201ll ) can easily be rewritten to show the uncertainty factor U' remaining 



in the rate after a measurement, 

U' = U!^^X + {1-X)K = U!^^ + {1-X){U-U!^^) , (12) 

where iYfxp is the experimental uncertainty factor. A perfect measurement without any error bar 
would have UHp = 1, whereas a 30% uncertainty, for example, would t ranslate into UHp = 1.3. 



The uncertainty in X itself has been discussed in iRauscher et al.l (|201ll ). It scales inversely with 
the value of X and has the property to keep large values of X large and small ones small, thereby 
leaving unchanged any conclusions regarding the impact of a measurement. 

Figures [2H7] show X for a number of reactions. Values for other reactions and temperatures 
can be extracted from Tables [U [21 and [3l As can be seen, the values depend on temperatur e and 



NLD for a given reaction, not on separation from the line of stability. iRauscher et al.l (|201ll ) have 
pointed out that even neutron captures on intermediate and heavy stable nuclei exhibit small X 
at s-process temperatures and in consequence many s-process reactions are less well constrained by 
experiments than often assumed. 

If X « 1, a measurement of the cross section would be the preferred way to proceed but 
it has to be performed in t he energy r ange required for the calculation of the astrophysical rate. 



These energies are given in IRauscher! (j20ld ). If it is not possible to measure at these energies, 
experiments can only indirectly - with the involvement of theory - improve a rate. If a measurement 
at higher energies is possible. Table H] has to be consulted to check up to which energy the cross 
section is sensitive to the same widths as the rate. If the accessible cross sections exhibit different 
dependences, a cross section measurement is futile and alternative approaches have to be considered 
to improve the prediction of a rate. The alternative is to provide data to improve the calculation 
of the actually relevant widths. Table [7] shows which input is mainly relevant for which type of 
widths. Measurements can be performed to either directly improve the input or to provide data to 
test and improve the models used to calculate the input. Experimental determination of masses and 
nuclear spectroscopy for low- lying states are examples for the direct improvement, measurements of 
scattering and reaction cross section for a better determination of optical potentials are examples 
for providing local and global data for improving input to the calculation of the rate. 

By examining Tables [IHSl theorists can directly see which nuclear properties affect an astro- 
physical rate most and may be in need of improvement. They can also immediately see the impact 
a different model prediction has on reaction rates and cross sections. One issue, already pointed out 
in §[11 should be remembered in this context. A similar procedure as suggested here for judging the 
sensitivities and resulting uncertainties in the rates has also to be applied to each input stemming 
from a theoretical model. The sensitivities of that model to its input have to be determined sep- 
arately to be able to arrive at a systematic evaluation of the connected uncertainties. Comparing 
the results of a limited selection of models is not a systematic variation and might lead to incorrect 
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conclusions when models fortuitously obtain similar values for astrophysically relevant quantities. 



2007; 



astrophysical reaction rates (see, e.g., 


Gviirkv et al. 


2006, 


2010: 


Rauscher 


2011 
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Sauerwein et al. 


2011: 


Halasz et al. 



20121 ). Such plots can be made using the data in Tables [IHSl they are also available online at 



|http : //nucastro . org/nonsmoker . html, Sample codes to extract plotting data for single reac- 
tions from the tables included here can be found at http : / / download . nucastro . org/codes/ . 



5. Summary 

Sensitivities of astrophysical reaction rates to systematic variations of reaction widths and the 
NLD have been presented and general dependences have been discussed. Additionally, the sensitiv- 
ities of reaction cross sections to the same quantities have been investigated. Considered were the 
compound nucleus reactions (n,7), (n,p), (n,a), (p,7), (p,n), (p,a), (0,7), (a,n), (a,p), (7,n), (7,p), 
and (7,0) for target nuclei between the driplines and 10 < Za < 83. General sensitivity trends 
have been discussed in detail and recommendations for application of the sensitivity information 
to assess and constrain uncertainties in astrophysical reaction rates were given. 

The sensitivities and the shown trends are robust results which do not largely depend on the 
specific reaction code or the input used because they are calculated from ratios of quantities (rates 
or cross sections) and the basic dependence of the result on a given property is determined by funda- 
mental nuclear properties and should be independent of the model code used. The only restriction 
is that a compound reaction mechanism is dominating the rate, either proceeding through narrow, 
discrete resonances or a large number of unresolved resonances. The (simpler) determination of 
sensitivities in the direct reaction mechanism has been outlined in § [3l 

This research was supported in part by the European Commission within the FP7 ENSAR/THEXO 
project and by the EuroGENESIS Collaborative Research Programme. 
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Fig. 1. — Energy scheme of a compound reaction a + A ^ C ^ B + b, proceeding via a compound 
state at energy Eq with spin Jq and parity ttc- Transitions to and from states in the target 
nucleus and Bi, in the final nucleus arc shown. In a capture reaction, the "ejectile" is a photon 
and nuclei B and C coincide. In a laboratory reaction, only the transition from the g.s. (oq) 
contributes in A whereas all transitions to and from the g.s. and all excited states of A and B with 
Efx,u < Be — Spro,ei, respectively, contribute in the effective cross section a^^ . 
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Fig. 2. — Ground state contribution X to stellar (n,7) rates at 1.5 GK. See the electronic edition 
of the Journal for a color version of this figure. 



- 25 - 



90 I 1 1 1 1 1 1 1 1 1 1 I — I 1 




20 40 60 80 100 120 140 160 180 200 
neutron number 



Fig. 3. — Ground state contribution X to stellar (7,n) rates at 1.5 GK. See the electronic edition 
of the Journal for a color version of this figure. 
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Fig. 4. — Ground state contribution X to stellar (p,7) rates at 1.5 GK. See the electronic edition 
of the Journal for a color version of this figure. 
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Fig. 5. — Ground state contribution X to stellar (7,p) rates at 1.5 GK. See the electronic edition 
of the Journal for a color version of this figure. 
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Fig. 6. — Ground state contribution X to stellar (0,7) rates at 1.5 GK. See the electronic edition 
of the Journal for a color version of this figure. 
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Fig. 7. — Ground state contribution X to stellar (7,0) rates at 1.5 GK. See the electronic edition 
of the Journal for a color version of this figure. 
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Fig. 8. — Sensitivity of the stellar rate of ^^^Sm(a,7)^^^Gd to a variation of various widths as 
function of plasma temperature. 
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Fig. 9. — Sensitivity of the stellar rate of ^^^Sn(n,7)^"^^Sn to a variation of various widths as 
function of plasma temperature. There is zero sensitivity to both proton and a widths and their 
lines overlap. 
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Fig. 10. — Sensitivity of the stellar rate of ^^Ru(p,7)^''Rh to a variation of various widths as function 
of plasma temperature. The sensitivities to changes in both neutron and a widths are negligible 
and their lines overlap. See the electronic edition of the Journal for a color version of this figure. 
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Fig. 11. — Sensitivities of stellar proton capture rates at 1.5 GK to a variation of the proton width. 
See the electronic edition of the Journal for a color version of this figure. 
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Fig. 12. — Sensitivities of stellar proton capture rates at 1.5 GK to a variation of the 7- width. See 
the electronic edition of the Journal for a color version of this figure. 
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Fig. 13. — Sensitivities of stellar a capture rates at 1.5 GK to a variation of the 7- width. See the 
electronic edition of the Journal for a color version of this figure. 
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Fig. 14. — Sensitivities of stellar neutron capture rates at 0.3 GK to a variation of the neutron 
width. See the electronic edition of the Journal for a color version of this figure. 
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Fig. 15. — Sensitivities of stellar neutron capture rates at 0.3 GK to a variation of the 7-width. 
See the electronic edition of the Journal for a color version of this figure. 
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Fig. 16. — Sensitivities of stellar neutron capture rates at 0.3 GK to a variation of the proton 
width. See the electronic edition of the Journal for a color version of this figure. 
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Fig. 18. — Sensitivity of the laboratory cross section of ^^Ge(n,p)^^Ga as function of cm. energy. 
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Fig. 19. — Sensitivity of the laboratory cross section of ^^'^Hf(n,7)^^''"Hf as function of cm. energy; 
note the logarithmic energy scale. The astrophysically relevant energy range is 0. < £^ < 0.1 MeV 
for s-process nucleosynthesis. 
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Fig. 20. — Sensitivity of the laboratory cross section of ^^^Hf(7,n)^^°Hf as function of cm. energy. 
The sensitivities to variations of proton or a width are zero and not shown here. 
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Fig. 21. — Energy -Eiast of the highest included discrete, experimental level; above that energy, a 
theoretical NLD is used. See the electronic edition of the Journal for a color version of this figure. 
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Fig. 22. — Sensitivity of the stellar rate of ^^Ru(p,7)^^Rh to a variation of the NLD (with and 
without inclusion of experimental excited states) in the proton and 7 channel as function of plasma 
temperature. See the electronic edition of the Journal for a color version of this figure. 
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Fig. 23. — Sensitivity of the laboratory cross section of ^^Ru(p,7)^^Rh to a variation of the NLD 
(with and without inclusion of experimental excited states) in the shown channels. See the electronic 
edition of the Journal for a color version of this figure. 
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Fig. 24. — Sensitivity of the laboratory cross section of ^^^Sm(Q;,7)^^^Gd to a variation of the NLD 
(with and without inclusion of experimental excited states) in the shown channels. 
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Fig. 25. — Sensitivity of the laboratory cross section of i'^°Hf(n,7)^^^Hf to a variation of the NLD 
(with and without inclusion of experimental excited states) in the shown channels. See the electronic 
edition of the Journal for a color version of this figure. 



Table 1. Sensitivities of the astrophysical reaction rates to variations of different widths at 24 plasma temperatures; ground 

state contributions to the stellar rate are also given (see text for details) 



Element 




Ma 


Proj. 


Ejec. 










X 










X 














T 


= 0.1 GK 








T = 10.0 GK 




mo 


42 


92 


n 


g 


0.56 


0.35 


0.00 


0.00 


1.000 


0.87 


0.07 


-0.02 


0.00 


1.02 X 10-1 


mo 


42 


92 


n 


P 


-0.37 


0.47 


1.00 


0.00 


1.000 


-0.01 


0.08 


0.79 


-0.01 


6.57 X 10"'' 


mo 


42 


92 


n 


a 


-0.25 


0.31 


0.00 


1.00 


1.000 


0.00 


0.11 


-0.15 


0.96 


6.09 X 10-5 


mo 


42 


92 


P 


g 


0.00 


0.00 


1.00 


0.00 


1.000 


1.00 


-0.06 


0.09 


0.00 


1.56 X 10-1 


mo 


42 


92 


a 


g 


0.00 


0.00 


0.00 


1.00 


1.000 


0.65 


-0.49 


-0.10 


0.96 


3.86 X 10-1 



Note. — Table [T] is published in its entirety in the electronic edition of the Astrophysical Journal Supplement. A portion is shown here for 
guidance regarding its form and content. 
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Table 2. Ground state contributions X to the stellar rate at 24 plasma temperatures for 
reactions with negative Q value (except captures; see text for details) 



Element 


Za 


Ma 


Proj. 


Ejec. 


X 


X 


X 


X 












T = 0.1 GK 


T = 0.15 GK 


T = 9.0 GK 


T = 10.0 GK 


mo 


42 


92 


P 


n 


3.73 X lO-** 


3.73x10-'' 


2.53 X 10-3 


1.39 X 10-3 


mo 


42 


92 


P 


a 


1.00 X 10° 


1.00 X 10° 


2.98 X 10"* 


1.09 X 10-4 


mo 


42 


92 


a 


n 


1.00 X 10° 


1.00 X 10° 


6.27 x 10-2 


2.63 X 10-2 


mo 


42 


92 


a 


P 


1.00 X 10° 


9.98x10-1 


3.25 x 10-2 


1.17 X 10-2 


mo 


42 


93 


P 


n 


8.72 X IQ-^ 


8.84x10-1 


3.36 X 10-2 


1.53 X 10-2 



Note. — Table [2] is published in its entirety in the electronic edition of the Astrophysical Journal Supplement. A 
portion is shown here for guidance regarding its form and content. 



Table 3. Ground state contributions X to the stellar rate at 24 plasma temperatures for 
photodisintegration reactions (see text for details) 



Element 


Za 


Ma 


Ejec. 


X 


X 


X 


X 










T = 0.1 GK 


T = 0.15 GK 


T = 9.0 GK 


T = 10.0 GK 


mo 


42 


92 


n 


5.93 X 10-^ 


5.93x10-^ 


4.14 X 10-* 


3.32 X 10-* 


mo 


42 


92 


P 


2.03 X 10-* 


2.03x10-* 


2.09 X 10-3 


1.56 X 10-3 


mo 


42 


92 


a 


1.00 X 10" 


1.43x10-1 


1.48 X 10-3 


6.36 X 10-* 


mo 


42 


93 


n 


3.92 X 10-2 


3.92x10-2 


3.13 X 10-3 


2.14 X 10-3 


mo 


42 


93 


P 


2.15 X 10-^ 


2.15x10-^ 


3.50 X 10-* 


2.69 X 10-* 



Note. — Table [3] is published in its entirety in the electronic edition of the Astrophysical Journal Supple- 
ment. A portion is shown here for guidance regarding its form and content. 



Table 4. Sensitivities of laboratory cross sections to variations of different widths (see text for details) 



Element Za Ma Proj. Ejec. riEn E (MeV) """s^ '^^Sn '^^Sp ==Sa • ■ • E (MeV) ^^^Sg ^^^Sn ''"sp '^^Sa 



mo 


42 


92 


n 


g 


50 


0.001 


0.44 


0.49 


0.00 


0.00 


14.131 


0.99 


0.15 


-0.22 


-0.03 


mo 


42 


92 


n 


P 


50 


0.001 


-0.63 


0.74 


1.00 


0.00 


14.131 


0.00 


0.15 


0.64 


-0.03 


mo 


42 


92 


n 


a 


50 


0.001 


-0.30 


0.37 


0.00 


1.00 


14.131 


0.00 


0.15 


-0.22 


0.95 


mo 


42 


92 


P 


g 


50 


0.270 


0.00 


0.00 


1.00 


0.00 


7.141 


1.06 


0.00 


0.00 


0.00 


mo 


42 


92 


P 


n 


50 


8.626 


0.00 


0.00 


0.00 


0.00 


14.141 


-0.01 


0.21 


0.52 


0.00 



Note. — Table |4] is published in its entirety in the electronic edition of the Astrophysical Journal Supplement. A portion is shown here for guidance 
regarding its form and content. 



Table 5. Sensitivities of the astrophysical reaction rates to variations of nuclear level densities (with and without the 
inclusion of discrete excited states) in the 7-, neutron-, proton-, and a-channel at 24 plasma temperatures (see text for details) 



Element 


Za 


Ma 


Proj. 


Ejec. 


^last 




last 


-^last 


13 




13 *P 




r „0 
D*g 


r „0 


D*P 


D*a 


























T = 


0.1 GK 








mo 


42 


92 


n 


g 


2.77 


5.12 


0.501 


4.10 


0.46 


0.00 


0.00 


0.00 


0.55 


0.00 


0.00 


0.00 . . . 


mo 


42 


92 


n 


P 


2.77 


5.12 


0.501 


4.10 


-0.31 


0.00 


0.00 


0.00 


-0.11 


0.00 


1.00 


0.00 . . . 


mo 


42 


92 


n 


a 


2.77 


5.12 


0.501 


4.10 


-0.20 


0.00 


0.00 


0.00 


-0.24 


0.00 


0.00 


0.01 ... 


mo 


42 


92 


P 


g 


1.79 


0.966 


5.12 


0.659 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 . . . 


mo 


42 


92 


a 


g 


3.95 


1.35 


2.23 


5.12 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 . . . 



Note. — Table [5] is published in its entirety in the electronic edition of the Astrophysical Journal Supplement. A portion is shown here for guidance 
regarding its form and content. 



Table 6. Sensitivities of laboratory cross sections to variations of nuclear level densities in the 7-, neutron-, proton-, and 
a-channel, with and without inclusion of discrete excited states (see text for details) 



Element Za Ma Proj. Ejec. riEn E (MeV) gsg gsn gsp g^a g^O gs^ gsg 



mo 


42 


92 


mo 


42 


92 


mo 


42 


92 


mo 


42 


92 


mo 


42 


92 



n g 50 

n p 50 

n a 50 

P g 50 

p n 50 



0.001 


0.43 


0.00 


0.001 


-0.45 


0.00 


0.001 


-0.22 


0.00 


1.378 


0.00 


0.00 


8.626 


0.00 


0.00 



0.00 


0.00 


0.51 


0.00 


0.00 


-0.21 


0.00 


0.00 


-0.27 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 



0.00 0.00 0.00 

0.00 1.00 0.00 

0.00 0.00 0.01 

0.00 0.00 0.00 

0.00 0.00 0.00 



Note. — Table [6] is published in its entirety in the electronic edition of the Astrophysical Journal Supplement. A portion is shown here 
for guidance regarding its form and content. 
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Table 7. Input to the compound reaction calculation and its impact (in no particular order). 

Quantity Impact 

Nuclear masses all widths'^ 

Properties of g.s. and excited states particle widths'' 

Nuclear level density 7-width'' 

Optical potentials particle widths'^ 

7 strength function 7-widths'^ 

Deformations optical potentials, NLD, 7 strength function 

(Nuclear matter density distribution) (optical potentials)*^ 

'^Determining the relative energy in all transitions through separation energies 

''In principle, all widths are affected by inclusion of discrete levels up to an energy 
-^last) with specific excitation energy, spin, parity, and the use of a NLD above, see 



Equation (lllh . Relevant 7-energies, however, are 2 — 4 MeV and mainly lead to states 
above -Eiast whereas particle transitions prefer low-lying states. 

^Separate optical potentials per particle type 

'^El is dominating, depending on NLD sometimes also Ml relevant 

'^Some optical potential descriptions require the knowledge of nuclear densities, e.g., 
foldir ig potentials or microscopic potentials based on the local density app roximation 



(e.g., iJeukenne. Leieune. &: Mahauxlll977l : iBauge. Delaroche. &: Girodll200ll ) 



